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Abstract 

Using projector augmented wave approach within the generahzed gradient approximation, we 
have studied the structural property and electronic structure of ZrBi2. The calculated lattice 
constants and bulk modulus are in good agreement with the available experimental values. A 
detailed study of the electronic structure and the charge-density redistribution reveals the features 
of strong covalent B-B and weak covalent Zr-B bondings in ZrBi2. The states at the Fermi level 
mainly come from the cr-Py and vr-p^ orbitals of boron atoms, which are slightly hybridized with 
the t2g-dxz and t2g-dyz orbitals of Zr atoms. As the increased hydrostatic pressure on ZrBi2, the 
total density of states at the Fermi level decreases. 
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I. INTRODUCTION 



Metal borides have been the subject of extensive research in recent years |T], ^ |^ because 
they exhibit unique properties such as high mehing points, high hardness values, and excel- 
lent oxidation resistance. These particular properties mean borides are promising materials 
for new heat-resistant, corrosion-resistant, and wear-resistant alloys and coatings 0]. In par- 
ticular, the discovery of superconductivity (Tc ^ 39 K) in magnesium diboride (MgB2) 0, 
which is crystallized in a simple layered hexagonal AIB2 structure, has stimulated a re- 
newed interest to search for new superconductors and to understand their mechanism of 
superconductivity in borides. 

ZrB 12 is found to have the highest superconductivity critial temperture {Tc ~ 6 K) in 
the family of dodecaborides MB12 H, 0, In order to understand its superconductivity, 
some efforts are made by measuring the specific heat, magnetic suseptibility, resistivity, 
thermal expansion, electron transport, penetration depth, and the upper critical magnetic 
field of ZrBi2 0, ||, |, 0, 0, [l^, 0, |16|- Specially, a negative pressure effect on 



Tc of ZrBi2 is observed in its magnetization measurements under pressure, which contrasts 
to the previous suggestion that the superconducting critical temperature of ZrBi2 might 



increase under pressure [17]. The previous band structure calculations on ZrBi2 mainly 



concern the electronic structure of ZrBi2 under ambient pressure [0, |I9|. Therefore, we 
perform ab initio total energy calculations to study the properties of ZrBi2, epescially about 
its electronic structure under high pressure. 

II. METHOD 

All calculations on ZrBi2 are performed using the projector augmented wave (PAW) 



approach as implemented in the VASP code [20, 21 1. The exchange-correlation energy 



functional is treated with the generalized gradient approximation (GGA) with Perdew and 



Wang parameterization (known as GGA-PW91) [22|. Electron-ion interaction is represented 



by the PAW method [^, ^ and wave functions are expanded by the plane waves up to 
an energy cutoff of 500 eV. Brillouin-zone integrations are approximated using the special 
A;-point sampling of Monhkorst-Pack scheme [Q with a 9x9x9 grid. The above calculation 
parameters of energy cutoff and /c-mesh make the total energy convergence to 2 meV/atom. 
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The internal atomic coordinates of cubic ZrBi2 are relaxed at a series of fixed volumes until 
the force exerting on each atom is less than 0.01 eV/A. The obtained energies are fitted 
with the third-order Birch- Murnaghan equation of state (EOS) to give the equilibrium 



volume and the minimum energy. The final calculated cell parameters are listed in Table |. 

III. RESULTS AND DISCUSSIONS 

The crystalline structure of ZrBi2 is cubic and belongs to space group Of^-FmSm (No. 225) 
p7| , shown in Fig. 0. The Wyckoff positions of Zr and boron atoms are 4a (0, 0, 0) and 
48z (1/2, X, x), respectively, i.e. the Zr atoms and cuboctahedral B12 clusters are arranged 
as in the NaCl-like structure. The 13-atom fee unit cell of cubic ZrBi2 is characterized 
by single lattice constant and the internal parameter of boron atoms. In present work, 
we obtain the lattice constant a is 7.409 A and the internal coordinate parameter x of 
boron atoms is 0.161, in good agreement with experiment values |2^ (See Table |). 



Based on GGA calculations, the bulk modulus and its pressure derivative for ZrBi2 are 
predicted as 233.0 GPa and 3.64, respectively. Note that in the investigation of electron- 
phonon coupling depentent on pressure of ZrBi2 the authors used the bulk modulus of an 
analogy compound UB12 {Bq = 249 GPa) because the experimental value of bulk modulus 
of ZrBi2 is not reported. In addition, the predicted bulk modulus of ZrBi2 is close to that of 
ZrB2 (245 GPa obtained from the ultrasonic measurements p9| or 238 GPa predicted from 



ab initio calculations ||3^) crystallized in hexagonal layered AlB2-like structure. 

The energy band structure and electronic density of states (DOS) of cubic ZrBi2 in 
equilibrium are shown in Fig. ^(a) and ^, respectively. These may be compared well with 



previously calculated results of ZrBi2 based on a variety of computational methods pq , |19 
Overall these suggest that ZrBi2 exhibits a metallic character with the mixed p-d states 
rather than a simple p- or metal [^J. The valence bands in Fig. ^(a) are split into three 
disjointed groups. The lowest group centered at -15 eV below the Fermi level (Ep) has a 
width of about 1 eV and the corresponding bands mostly come from the bonding states 
of B-2s orbitals. Combining from the DOS figures (^), the higher one in the energy range 
from -13 eV to -10 eV below the Ep is largely contributed by the antibonding states of B-2s 
orbitals. The upper one in the energy range from -10 eV to Ep is mainly composed by 
B-2p orbitals and Zr-Ad electrons. In addition, it suggests that there is a weak hybridization 
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between B-2p and Zi-Ad orbitals occuring below the Ep. The total density of states at the 
Fermi level N{Ef) is about 1.478 states/eV-unit cell, in consistence with the calculated value 
in other work |T^. The states at Ep is mainly made by the 2p electrons of boron atoms 
{cr-py and Tc-p^ orbitals) and 4d levels of Zr atom {t2g-dyz and t2g-dxz orbitals). This can 
also be seen from the the partial charge density calculated within a 0.2 eV energy window 
around Ep for ZrBi2, which is shown in Fig. ^(b). 

In order to understand the Zr-B and B-B bondings of ZrBi2 in equilibrium, the line 
charge densitiy along nearest neighbor Zr— B and B— B pairs are illustrated in Fig. ^. The 
calculated bond lengths of atoms are dzi-B = 2.75 A and c^b-b = 1-68 A, which agree well 



with the experimental values It is clearly seen that charge highly accumulates around 
Zr atom due to its more valence electrons than that of boron atom. Contrast to the strong 
covalent character of B-B bonding, the Zr-B bonding exhibits a weak covalent character. It 
also supported by above discussion that the d orbital of Zr atoms and p orbital of boron 
atoms have a weak hybridization. In order to further reveal the topology of the Zr-B and 
B-B bondings. Fig. ^(b) illustrates the contour plots of electron density on crystallographic 
(010) plane. It again indicates that B-B bonding exhibits a strong covalent character. 

Now we turn to discuss the effect of high hydrostatic pressure on electronic structure of 
ZrBi2. The values of DOS at the Ep for ZrBi2 under a serial of hydrostatic pressure are 
listed in Table y. It can be seen that as the pressure increased from to 22.1 GPa, the 
total density of states at the Ep decreases about 0.264 states/eV-unit cell. Considering the 
electron-phonon coupling constant Aei-ph depends linearly on the N{Ep) [^, our calculated 
results indicate that the Aei-ph of ZrBi2 would decrease under the high pressure. It also 
supported by recent experimental observation that there is a negative pressure effect 
on Tc and Aei-ph of ZrBi2 with dTc/dp of -0.0225 K/kbar and dlnX^i-ph/dp of -0.2%/kbar. 
Additionaly, as the hydrostatic pressure increased, the contribution to N{Ep) of Zr atoms 
decrease than that of boron atoms. This could be understood from that the B-B bonding is 
stronger than Zr-B bondingbond, and hence the bond length of B-B atoms decreases more 
slightly under the high pressure (see Table ||). Upon the hydrostatic pressure to 22.1 GPa, 
the width of valence bands of ZrBi2 is increased to -16 eV as shown in Fig. |^(b), companied 
with the 2s-bonding states of B boron atoms toward lower energy. 
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IV. CONCLUSIONS 



The structural property and electronic structure of ZrBi2 have been calculated and an- 
alyzed. Our calculated structural parameters including the lattice constants and internal 
parameters are in good agreement with the experiment data. The detailed study of the 
electronic structure and charge density redistribution reveals the strong covalent of bonding 
of B-B atoms and role of p orbitals of boron atoms and d orbitals of Zr atom in the super- 
condutivity of ZrBi2. Our calculated results also indicate that there is a negative pressure 
effect on the superconductivity of ZrBi2. 
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TABLE I: Calculated lattice parameter (a, in A), atomic coordinate parameter {x for boron atoms), 
volumes of unit cell (Vq, in A^), and bulk modulus {Bq, in GPa) of cubic ZrBi2, along a comparison 
with other theoretical work and available experiment values. 



Bn 



This work 
Other work IT 



Other work |1S] 



Expt. 123] 
Expt. pi 



7.409 
7.419 
7.408 
7.408 
7.408 



0.1696 

0.1696 
0.166 
0.1699 



101.69 
102.09 
102.09 
101.63 
101.63 



233.0 



TABLE IL The bond lengths of atoms {dzr-B and d^-B, in A), total and partial density of states at 
the Fermi level {N{Ef) and NpariEp), in states/eV-unit cell) of cubic ZrBi2 at different hydrostatic 
pressures. 



Pressure 


dzi-B 


^B-B 


N{Ef) 


NpariEp) 

Zr-4d 


B-2p 


GPa 


2.75 


1.68 


1.478 


0.360 


0.585 


4.9 GPa 


2.73 


1.67 


1.306 


0.315 


0.518 


10.2 GPa 


2.71 


1.66 


1.251 


0.298 


0.503 


22.1 GPa 


2.68 


1.64 


1.214 


0.284 


0.501 



Figure Captions 

• FIG |l|: The ball-and-stick model for unit cell of ZrBi2. Zr atoms are shown as dark 
gray balls and boron atoms as light gray balls. 

• FIG 0: Electronic band structure of ZrBi2 (a) in equilibrium and (b) at 22.1 GPa 
hydrostatic pressure. The zero of energy is set as the Fermi level and shown in dash 
line. 

• FIG 1^: (a) Total and partial density of states (DOS) , (b) and (c) angular momentum 
decompsed DOS of ZrBi2 in equilibrium. 

• FIG ^: Contour plots of the charge density on (010) surface of ZrBi2 in equilibrium: 
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FIG. 1: 



(a) partial charge density with the energy window of [Ep — 0.1, Ep + 0.1] is in an 
increment of 2.0x10"^ e/A^ from to 2xl0~^ e/A^ and (b) total charge density is in 
an increment of 0.1 e/A^ from to 1.5 e/A^. 

• FIG H: Charge density at the equilibrium lattice parameter of ZrBi2, along the lines 
strarting from: (a) Zr atom to nearst neighbor B atom, and (b) B atom to nearest 
neighbor B atom. 
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